Resonance decay effects on anisotropy parameters 
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Abstract 

We present the elliptic flow V2 of pions produced from resonance decays. The transverse mo- 
mentum px spectra of the parent particles are taken from thermal model fits and their V2 are 
fit under the assumption that they follow a number-of-constituent-quark (NCQ) scaling law ex- 
pected from quark-coalescence models. The Vi of pions from resonance particle decays is found to 
be similar to the measured pion V2- We also propose the measurement of electron V2 as a means 
to extract open-charm V2 and investigate whether a thermalized system of quasi-free quarks and 
gluons (a quark-gluon plasma) is created in collisions of Au nuclei at RHIC. 
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Introduction. — One of the surprising obser- 
vations made at RHIC is the measurement 
of a number-of-constituent-quark (NCQ) 
dependence for both elliptic flow v 2 and 
the nuclear modification factor Rcp at in- 
termediate pt (1.5 < pr < 5 GeV/c) [1]. 
Models of hadron formation by constituent- 
quark coalescence provide a viable expla- 
nation for these observations whereas ex- 
pectations based on conventional fragmen- 
tation approaches are inconsistent with 
the data [2,3,4]. In coalescence models an 
NCQ-scaling of t>2 arises as a consequence 
of hadrons coalescing out of a thermal 
distribution of partons and reveals the 
flow developed during a partonic epoch at 
RHIC. Pion v 2 , however, appears to violate 
NCQ-scaling. In this paper, we study the 
effect of resonance decays on pion v 2 - We 
show that when decays are taken into ac- 
count, the measured pion i; 2 may become 
consistent with the NCQ-scaling demon- 



strated by kaon (K + , K , Kg), proton, A, 
and H v% [1,5]. 

The particle azimuthal distribution with 
respect to the reaction plane at rapidity y 
can be described by a Fourier expansion: 



dN 

~dK4> 



oc 1 + S„2t> n cos(nA0), 



where is the difference in azimuth an- 
gle between the particle and the reaction 
plane. The first and second Fourier coeffi- 
cients, v\ and V2, historically are called di- 
rected and elliptic flow, respectively. All co- 
efficients can be calculated from the rela- 
tion: v n = (cos(nA0)). 

As the volume of the system created in 
an off-axis nucleus-nucleus collision ex- 
pands, its spatial anisotropy quenches. 
The momentum-space anisotropics repre- 
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sented by the Fourier coefficients v n pre- 
serve information about the early collision 
dynamics when the spatial anisotropy was 
largest [2,6,7,8]. Since the initial overlap 
region is elliptical in shape, the second har- 
monic coefficient v 2 is the largest and most 
studied. 

NCQ-scaling ofv 2 . — Fig. 1 shows the 7r + + 
7r~, Kg, p+p, and A + A v 2 from minimum- 
bias 197 Au + 197 Au collisions at ^/s^n = 
200 GeV [1,5]. In the lower pr region (px < 
1.0 GeV/c), the values of v 2 are lower for 
higher mass hadrons. Hydrodynamic calcu- 
lations [9] predict the observed mass depen- 
dence of v 2-perhaps implying that a ther- 
malized system has been created in colli- 
sions at RHIC energy. At higher p? {pr > 
2 GeV/c), the v 2 measurements saturate at 
values below the hydrodynamic model pre- 
dictions. The saturated value of v 2 and the 
Pt scale where the saturation sets in de- 
pends on the particle-type: the baryon v 2 
saturates at higher pt and at larger values 
than that of mesons. 

According to coalescence models [4], after 
scaling both v 2 and pt with the number of 
the constituent quarks (NCQ) in the cor- 
responding hadron, all particles at inter- 
mediate pt should fall onto one universal 
curve. The NCQ-scaled v 2 measurements 
in Fig. l-(b) show that v 2 (pt /n q ) /n q for 
PT/n q > 0.6 GeV/c is similar for all parti- 
cles except pions. This observation, coupled 
with the NCQ-dependence observed at in- 
termediate pt in the nuclear modification 
factor Rcp is evidence of hadron formation 
by coalescence or recombination. In this 
case, v 2 (pT/n q )/n q represents a constituent 
quark momentum-space anisotropy v\ that 
arises as a consequence of collectivity in a 
partonic stage. Based on coalescence mod- 
els, NCQ-scaling suggests the creation of a 
quark-gluon plasma (QGP) with v\ char- 
acterizing the properties of the QGP. For 
this reason, understanding the source of the 
discrepancy in the NCQ-scaled pion v 2 is 
imperative. 



that NCQ scaling is valid for all hadrons 
other than pions and use the published 
v 2 measurements [1,5] to parameterize 
v 2 (pt l n q) / n q- The pt distributions are as- 
sumed to follow an exponential form with 
slope parameters taken from measurements 
when available. We use chemical fits to 
fix the relative hadron abundances [10,11]. 
Since the pion mass is much smaller than 
the sum of its constituent quarks masses, 
direct pions are not necessarily assumed 
to follow the scaling predicted from coales- 
cence models. As such, we do not consider 
direct pions, and instead choose to study 
the t>2 of the secondary pions. Given the 
model uncertainties, extraction of the di- 
rect pion v 2 is difficult and remains an open 
question. Finally we will discuss how to ex- 
tract open-charm v 2 based on the decayed 
electrons. 
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Fig. 1. (a) Measurements of the pt depen- 
dence of the event anisotropy parameters for 
7r, Kg, p, A. Dot-dashed lines are the results of 
fits; (b) Number-of-constituent-quark (NCQ) 
scaled v 2 . All particles except the pions follow 
the NCQ scaling. 



With this goal in mind, we study the ef- 
fect of secondary pions (from particle de- 
cays) on the measured pion v 2 . We assume 



Simulation results. — The v 2 values of the 
simulated resonances are parameterized by 



2 



fitting K s and A v-i [1] with the equation: 
1 + exp(— (pr/n — b)/c) 

(2) 

where a, 6, c and (i are the fit parameters, 
n is the constituent-quark number, and 
Pt is in the unit of GeV/c. The fit results 
are shown as dot-dashed lines in Fig. 1, 
where the fitting parameters are a = 0.1, 
b = 0.35, c = 0.2 and d = 0.03. The NCQ- 
scaling of v<i works well for kaons, protons 
and Lambdas within 0.6 < PT/ n q < 1-5 
GeV/c whereas pion w 2 deviates from NCQ 
scaling for all pt- The parameters from 
chemical fits are listed in Table 1. 

Table 1 

Parameters for the input resonances: the units 
for slope parameters T are GeV. The fraction 
of the hadrons are fixed from the measured 
abundances. 





Tl 


T2 


T3 


Fraction (%) 


p 


0.5 


0.4 


0.3 


60 ± 10 


U) 


0.5 


0.4 


0.3 


30 ± 10 




0.3 


0.3 


0.3 


6±5 


K* 


0.5 


0.4 


0.3 


2± 1 


A 


0.55 


0.55 


0.55 


2±2 



In high-energy collisions, a large fraction 
of hadrons are produced through resonance 
decays. This is particularly true for pions 
in high-energy heavy-ion collisions. At mid- 
rapidity, in collisions at RHIC, as many 
as 80% of pions are from resonance decays 
[12]. The dominant decays are p — > nn, 
uj -> 3tt, AT* (892) -> Ktt, K° s -> rnr and 
A — ► Ntt. With such a potentially large 
fraction of pions arising from decays, ac- 
counting for their effect on the observed 
pion V2 is very important. 

The pt distributions of pions from reso- 
nance decays are shown in Fig. 2. Many pi- 
ons at low pt are generated from resonances 
at larger pt where presumably the parent 
particle v 2 is larger. As a result the decayed 
pions take on a relatively large v% value. 
The decays from the p- and tu-mesons dom- 
inate the secondary pion pt spectrum. In 
this plot, a slope parameter of T = 0.4 GeV 



is used for the p distributions. In periph- 
eral collisions, the STAR measured slope 
parameter is 319±4(stat.)±32 (syst.) MeV 
[13]. The simulated results are in a good 
agreement with the PHENIX 7r° data from 
minimum bias 197 Au+ 197 Au collisions [14]. 
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Fig. 2. Resonance decayed pion distribu- 
tions. The summed spectrum is shown as a 
dashed-line. In this simulation, a slope param- 
eter of T = 0.4 GeV is used for p,u> and K*. 
For Kg and A the respective slope parame- 
ters T = 0.3 GeV and 0.55 GeV are used. The 
relative fraction of the hadrons are listed in 
Table 1. For comparison, PHENIX 7r° results 
are shown as open-crosses. 

In Fig. 3, the v 2 values for the simulated 
decay pions are shown as dashed-lines. The 
resonances included in this study are the 
p, u, K*, Kg and A. The decay p — > ixix 
with a 100% branching ratio dominates the 
production of secondary-pions. The simu- 
lated resonance particles are restricted to 
mid-rapidity \y\ < 0.5. Increasing the ra- 
pidity window does not change the results. 
The pion v% in the region pt < 1.5 GeV/c is 
sensitive to the shape of the ppT spectrum. 
Dashed-, dotted-, and dot-dashed-lines cor- 
respond to the v 2 results from the different 
slope parameters listed in Table 1. For the 
smaller slope parameter T = 300 MeV, the 
decayed pion v 2 is below the data, leaving 
room for other contributions [15]. 

The effects of resonance decays on proton 
and kaon v% were also investigated. As in 
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Ref. [15], they were found to increase the 
observed Vi above the assumed NCQ-scaled 
input v 2 but by a much smaller degree than 
that of pions. Since the v 2 of kaons is used 
in our fits, resonance decays can change our 
assumed input t>2- These changes, however, 
will be small. In this paper we concentrate 
on the pion t> 2 . 




Transverse momentum p T (GeV/c) 
Fig. 3. The measured pion V2 (symbols) is com- 
pared to the simulated V2 for pions from reso- 
nance decays (dashed lines). The assumed V2 
of mesons and baryons are represented by the 
solid and dot-dashed lines, respectively. 

D-meson t>2- — NCQ scaling suggests that 
hadrons at intermediate pt are formed from 
a thermal partonic phase created in heavy- 
ion collisions at RHIC. In this system, the 
high initial matter density gradient and 
copious interactions among partons leads 
to a collective motion. The large v 2 values 
measured for multi-strange hadrons also 
indicate that partonic collectivity develops 
in collisions at RHIC [16]. The demon- 
stration of collectivity is necessary but 
not sufficient to show that local thermal 
equilibrium is established. If collectivity is 
developed by much heavier charm quarks, 
however, it suggests that interactions be- 
tween charm quarks and u-, d-, or s-quarks 
must have been frequent enough and strong 
enough for the lighter quarks to have ther- 
malized. For this reason, the measurement 
of heavy flavor (open-charm) t> 2 can probe 
the degree to which the lighter u-, d-, and 
s-quarks thermalize. 

Since a large fraction of heavy-flavor 
hadrons decay through leptonic modes, 
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Fig. 4. The V2 of electrons from D-meson de- 
cays (top) and 7T° decays (bottom). An input 
D-meson v% from a coalescence model with 
zero charm quark V2 [17] is shown as a solid 
line (top). The input meson V2 curve taken 
from our fit to the NCQ-scaled Kg and A + A 
V2 is shown as a dot-dashed line in both plots. 

electron V2 measurements may offer a con- 
venient way to study open-charm t> 2 . In Fig. 
4- (top), the V2 for electrons from D-meson 
decays is shown. Dot-dashed lines repre- 
sent the v 2 of the parent meson taken from 
our fit to the NCQ-scaled K s and A + A 
t> 2 . The solid line represents a D-meson v% 
where D-mesons are assumed to coalesce 
from light quarks with v%' > and charm 
quarks with v% — [17] . We simulate the D- 
meson px distribution using the PYTHIA 
event generator [18]. The shape of the cal- 
culated spectrum is well represented by a 
power-law function. Approximately 50 M 
D-meson events are used in this simulation. 
For electron px > 1.5 GeV/c, and a sat- 
urated v 2 at higher px the v 2 of electrons 
from D-meson decays becomes similar to 
the parent D-meson t> 2 . The decay elec- 
tron t> 2 is found to be sensitive to changes 
in the assumed D-meson t> 2 . As such, the 
measurement of electron v 2 can distinguish 

between v\ ~ t> 2 ' d and v\ ~ 0. The degree 
of heavy flavor thermal equilibrium can be 
assessed by measuring electron t> 2 within 
1 < Pr(e) < 3 GeV/c, a region which cor- 
responds to 2 < Pt(D) ^ 5 GeV/c. The 
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authors of Ref. [19] come to a similar con- 
clusion but do not evaluate background 
contributions to the electron spectrum. 

Neutral pion decays are the dominant 
source of "background" electrons. While 
the two-photon decay process, 
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can be identified by its decay topology 
[20] , the pion Dalitz decay can only be sub- 
tracted statistically. In Fig. 4-(b) we show 
the t>2 of electrons from simulated pion 
Dalitz decays. The pion distribution can 
be obtained from measurements at RHIC 
[14]. For these simulations a 100% con- 
version probability is assumed. In STAR, 
however, the probability is closer to 5%. 
The decayed electrons predominantly have 
p T < 0.5 GeV/c [20]. 

Electrons from heavy flavor decays begin 
to dominate the electron spectrum above 
Pt ~ 3 GeV/ c. With knowledge of the pion 
yield, D-meson yield, the pion v 2 , and the 
electron v 2 it will be possible to extract the 
D-meson t> 2 . These measurements can be 
made by both the PHENIX and STAR col- 
laborations at RHIC. Direct photon w 2 can 
also be measured with this method. 

Summary. — We have studied the effect of 
resonance decays on the pion v 2 in Au+Au 
collisions at RHIC. When the v 2 values 
for resonances are assumed to follow NCQ 
scaling, the pions generated in their de- 
cays take on v 2 values similar to those 
measured at RHIC. The dominant source 
of secondary pions is p decays. We have 
shown that when decays are accounted 
for, the measured pion v 2 values may be- 
come consistent with the NCQ scaling law 
that suggests the development of partonic 
collectivity in collisions at RHIC. Model 
uncertainties, however, make it difficult to 
extract the direct pion v 2 . In addition, we 
propose the measurement of electron t> 2 
as a means to study open-charm v 2 and 
hence the degree of thermalization reached 
at RHIC. 
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